Evidence for entangled states of two coupled flux qubits 
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We have studied the low-frequency magnetic susceptibility of two inductively coupled flux qubits 
using the impedance measurement technique (IMT), through their influence on the resonant prop- 
erties of a weakly coupled high-quality tank circuit. In a single qubit, an IMT dip in the tank's 
current -voltage phase angle at the level anticrossing yields the amplitude of coherent flux tunnel- 
ing. For two qubits, the difference (IMT deficit) between the sum of single-qubit dips and the dip 
amplitude when both qubits are at degeneracy shows that the system is in a mixture of entangled 
states (a necessary condition for entanglement). The dependence on temperature and relative bias 
between the qubits allows one to determine all the parameters of the effective Hamiltonian and 
equilibrium density matrix, and confirms the formation of entangled eigenstates. 

PACS numbers: 85.25.Cp, 03.67.Lx, 03.67.-a, 85.25.Dq 



While its exact role in the power of quantum comput- 
ers can still be debated, there is no doubt that entan- 
glement is a key feature of quantum registers, absent in 
their classical counterparts |l| . From the point of view of 
scalability, implementing quantum registers as integrated 
circuits is attractive. However, demonstrating entangle- 
ment in macroscopic solid-state systems remains a daunt- 
ing task, even after the observation of quantum coherence 
in many types of superconductingqubits 0, 0, 0, 03, @ 
and also entanglement in charge |7j and current-biased 
Josephson-junction [8| qubits, using pulse and spectro- 
scopic techniques respectively. In this paper we investi- 
gate two coupled three-junction flux qubits 0, 0| using 
a different method, providing a simple criterion for eigen- 
state entanglement, by studying the latter's influence on 
the two-qubit magnetic susceptibility. 

The impedance measurement technique (IMT) [l"H re- 
lies on monitoring the current and voltage in a high- 
quality, low-frequency resonant tank circuit inductively 
coupled to the qubit. In various modifications, this 
method was applied to map the Joscphson potential pro- 
file an d to observe Landau-Zener transitions [HI ] 
and Rabi oscillations 0] in a three-junction flux qubit. 
Due to the tank's weak coupling to the qubit and high 
quality factor, decoherence times can be as high as 
2.5 /is Note that the IMT approach requires a 
small decoherence rate compared to only the tunnel- 
ing amplitudes (0.5~1 GHz), not the tank frequency 
(10~20 MHz). This paper generalizes to two qubits an 
IMT study of coherent tunneling ^5|. The tank is fed 
a small-amplitude ac signal; its effective impedance, and 
therefore the current-voltage phase angle, are sensitive to 
the qubit susceptibility \. For, say, one qubit at tempera- 
ture T — 0, x 9%E~, the curvature of the ground-state 



energy £L vs external flux $ [a different expression for 
X is used in Eq. (@J below]. Thus, the level anticrossing 
due to flux tunneling is revealed by a large peak in \x\- 

Our system of two flux qubits inductively coupled to 
each other and to the tank is shown in Fig. The 
qubits were fabricated out of aluminum by conventional 
shadow evaporation, nominally 1 /im apart, at the cen- 
tre of a niobium pickup coil. The area of each qubit 
was 80 ^m 2 , with self-inductance L a / b 39 pH; the 
two biggest junctions in each qubit had critical current 
« 400 nA and Coulomb energy e 2 /2Ch w 3.2 GHz (for 
the third junction, these values are 10~20% smaller and 
larger, respectively). The mutual inductance between the 
qubits Mai, — 2.7 pH was estimated numerically from 
the electron micrograph. The Nb coil, which together 
with an external capacitor of Ct ~ 470 pF forms a par- 
allel tank circuit with inductance Lt ~ 130 nH, reso- 
nant frequency lot/2h — 20.139 MHz and quality factor 
Qt = ojtRtCt = 1680 (at 10 mK; Rt is an effective 
tank resistance) |l6j | , was fabricated using e-beam lithog- 
raphy. The external magnetic flux through the qubits 
was created by the dc component of the current in the 
coil /del , and by the bias current 7dc2 through a wire close 
to one of the qubits. This allowed independent control 
of the bias in each qubit. 

The system of Fig. ^is described by the Hamiltonian 
H = Ho+HT+Hi n t+Hdi ss , where the two-qubit Hamilto- 
nian in the two-state approximation is expressed through 
Pauli matrices as ^3 



H 



-A a a^-A b a^+e a 4 a He b a^ + Ja^4 b \ (1) 



Ht is the tank Hamiltonian (a harmonic oscillator), the 
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FIG. 1: Two-qubit system coupled to a resonant tank circuit, 
(a) Schematic, (b) Micrograph. 



qubit-tank interaction is 



(2) 



and -ffdiss describes the standard weak coupling of the 
qubits to a dissipative bath 0. Here It is the cur- 
rent through Lt- The coefficients are A Q /£, = Ma/^ ^Ia/b, 
where Ma/^T is the qubit-tank mutual inductance and 
I a /i, is the magnitude of the persistent current in the 
corresponding qubit. The qubit biases are given by 
e a = I a $o(fx - 5 + /shift), e 6 = I b $ (f K - \ + i]f shift ), 
where / x = ^/^o accounts for the external flux $ tx 7 dcl 
created by the Nb coil in both qubits, while the parame- 
ters /shift cx Idc2 and rj = Mb w /M aw < 1 describe the bias 
difference between the qubits created by the additional 
wire. Here M aw (Mf, w ) are the mutual inductances be- 
tween the a (b) qubit and the additional dc wire (for our 
sample, M avl and M& w were calculated numerically, yield- 
ing r\ = 0.32). The qubit-qubit coupling J = M a bI a Ib is 
positive because the two qubits are in the same plane side 
by side, leading to antiferromagnetic coupling (according 
to the north-to-south attraction law). 

The effective inductance of the tank, and therefore its 
resonant frequency, depend on the state of the qubits. 
In the IMT method , this dependence is observed by 
continuously monitoring the phase angle between the 
average tank voltage and bias current /bias (t) fed into the 
tank (Fig.nj. For small enough /bias, tan 9 is determined 
by the real part of the two-qubit contribution to the tank 



susceptibility \' at wt 0, viz., 

Qt 



tan 6 



xVt), 



(3) 



where x( w ) has a Kubo-type linear-response expression 
through retarded Green's functions of the qubit opera- 
tors tri"^ . For weak damping H diss , the latter can be cal- 
culated with the equilibrium density matrix correspond- 
ing to H in Eq. Q • It can be generally assumed that the 
eigenvalues of H , ^ — 1,2,3,4, are non-degenerate 
and the eigenstates normalized, (v\n) = 5^ u . Taking into 
account the qubits' interaction with a dissipative envi- 
ronment (lit |2fJ , a standard calculation gives [2l| 



Pp. - Pu 



Tilo + E n 



(4) 



where p M = ^^^/E/ - ^' 1 ] is the population of 
energy level /z, in the energy denominators are deco- 
herence rates, and the real matrix elements are 

R^ = ^</*ki a) HH^ 0) lA*> + A?(MkW|i/)<i/kW|A*> 

+ A a A h <^V)Ha( fc V} 

+ \ a \M^ ) W){v\° i ? ) \p)- (5) 

Substitution into Eq. J2J yields 

tane = -2^y p » 
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At low frequencies ~ wt <C \Efj,—E v \/h and for weak 
damping T^ v <C \E^—E v \/h, the do not affect tan6, 
but are responsible for establishing the equilibrium dis- 
tribution. 

The first two terms in Eq. JSJ) are clearly positive, 
and non-zero even if the two-qubit states are factorized. 
The first (second) term corresponds to the contribution 
of qubit a (b), and peaks near that qubit 's degeneracy 
point. These contributions are practically independent 
of whether the qubits' degeneracy points coincide or not. 

The last two terms in Eq. JSJ) describe coherent flipping 
of both qubits, which is only possible for non-factorizable 
(entangled) eigenstates \fi), \v). For J > 0, these terms 
are found to be negative. Therefore the difference (IMT 
deficit) between the coinciding two-qubit IMT dip and 
the sum of two single-qubit dips provides a measure of 
how coherent the two-qubit dynamics is, that is, whether 
entangled eigenstates of Hq are formed |2l| . Eigenstate 
entanglement is a necessary condition for the equilibrium 
state to be entangled; of course, the eigenstate popula- 
tions play a role as well. However, it is already a suffi- 
cient condition for performing quantum gate operations 
(if decoherence times are long enough), since one can also 
initialize the system in a nonequilibrium state, say, by let- 
ting it relax at a point in parameter space for which there 
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is a large gap above the ground state, and subsequently 
returning to the working point 0, Q • 

The measurements on the system of Fig. fallowed to 
determine all parameters of Hq. The data are presented 
in Fig. (tan 9 vs flux bias) and Fig. 0] (T-dependence 
of dip amplitudes). 

First, consider idc2 = 27.3 (—32.7) liA, i.e., unequally 
biased qubits. In this case, whenever one qubit is near- 
degenerate, the other is strongly biased, putting it into 
a unique classical ground state (at the low end of our 
T-range) in which it behaves trivially. Hence, one es- 
sentially observes single-qubit properties. We match the 
corresponding IMT dips of Fig. [3 with the predictions of 
Eqs. (JSJ), JSJ for both the width and amplitude; see Fig. [3] 
This yields A a /h = 550 MHz, A b /h = 450 MHz, and the 
persistent currents I a ~ lb = I p = 320 nA, for a flux- 
to-energy bias conversion factor I p $o/h = 990 GHz [22j |. 
The width of the dips is practically T-independent, as ex- 
pected . The T-dependence of the amplitudes (Fig. 0] 
squares and triangles) agrees with these values for A Q / b . 
(The saturation below 35-40 mK 0] is likely due to 
a discrepancy between the mixing-chamber and sample 
temperatures.) 

We can determine the coupling constant in Eq. Q 
from the single-qubit measurements, as J/h = 
M ab I^/h = 410 MHz. This agrees well with the value 
J/h = 420 MHz obtained by fitting the T-dependence of 
the coincident IMT dip amplitude above the saturation 
temperature (Fig. 0] circles) with the theoretical curve. 

Comparison of the single-qubit and coincident IMT 
dips shows that the contribution to tan O due to eigen- 
state entanglement is significant. Indeed, the amplitude 
of the central dip in Fig. [2 at T = 50 mK is 1.12, com- 
pared to the sum 1.69 of the single-qubit dips. This 
means that the entangled terms [last two in Eq. (0)] are 
responsible for a contribution w —0.57. As the dashed 
line in Fig. [3] shows, dropping these terms would yield 
the coincident IMT dip in marked deviation from both 
the full theory and experiment. 

The IMT deficit also confirms the sign of the qubit- 
qubit interaction: in the case of ferromagnetic coupling, 
the coincident dip should be larger than the sum of the 
single-qubit ones. 

Plotting the ratio of the coincident dip to the sum 
of single-qubit ones (Fig. ^ inset), we see that it grows 
with T, approaching w 1. This is to be expected: ther- 
mal excitations tend to destroy coherent correlations be- 
tween the qubits (which will then behave as independent 
quantum systems with zero IMT deficit) — eventually, of 
course, destroying the single-qubit IMT dips as well. 

At 50 mK, T is comparable to the characteristic qubit 
energies (at the two-qubit degeneracy point, the top ex- 
cited state is ~ 100 mK above the ground state). The 
measurement time equals at least the tank saturation 
time Qt/^t and is longer in practice, 1~10 ms. Since 
this exceeds any conceivable qubit relaxation time ~ , 
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FIG. 2: Normalized tangent of the current -voltage angle in 
the tank vs external bias "1>. From the lower to upper curves, 
the temperature of the mixing chamber is 10, 50, 90, 160 mK. 
A relative flux bias /shift between the qubits is created by 
changing the current Jdc2 in the additional wire. The shifted 
curves correspond to J dc2 = 27.3 (-32.7) fj,A. Two IMT dips 
are then observed, showing tunneling in each qubit at the cor- 
responding degeneracy point / x 7^ ^. The central curve has 
7dc2 = —2.7 /iA. Both qubits are degenerate simultaneously 
at / x = I. One IMT dip is observed, with an amplitude about 
33% (at small T) less than the sum of two separate dips. 
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FIG. 3: Theoretical fit of the normalized tangent of the 
current-voltage angle in the tank vs external bias Solid 
curves, from lower to upper: T = 50, 90, 160 mK. The in- 
teraction energy between the qubits is 420 MHz. The dashed 
central peak corresponds to T = 50 mK, but with the entan- 
gled contributions in © omitted. Relative flux biases are the 
same as for the experimental curves in Fig. [2] 



the system has time to equilibrate. Indeed, the good 
quantitative agreement between experiment (Figs. [2 |3J 
and theory [Eq. in a wide range of T confirms that 
our system is described by the equilibrium density ma- 
trix with the Hamiltonian JQ). In other words, it is in a 
mixture of entangled two-qubit eigenstates. 

A measure of the entanglement of an arbitrary state is 
the concurrence C |23| , changing monotonically from for 
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FIG. 4: Temperature dependence of the IMT dip ampli- 
tudes. Squares: qubit a; triangles: qubit b; circles: coinci- 
dent dip. The solid and dashed curves are theoretical fits 
(the dashed one is for qubit a; the 6-curve is essentially iden- 
tical for T > 50 mK). Inset: relative size of the coincident dip 
tan C / (tan <d a + tan Qb)', the dotted line is the theoretical 
curve. 



non-entangled to 1 for maximally entangled states. For 
a normalized pure state a|00) + f3\01) +7|10) + <5|11), one 
has C = 2\a8—j3 r )\ \ in general, C depends both on the ex- 
istence of entangled eigenstates and on their population, 
and can be erased by temperature faster than quantum 
coherence in individual qubits |24j . Using the experimen- 
tally established parameters of the Hamiltonian QJ, the 
concurrences of its eigenstates at the two-qubit degener- 
acy point are immediately calculated: C\ = C4 = 0.39, 
C2 = C3 = 0.97. For the equilibrium density matrix of our 
system at 10 mK, these give C oq = 0.33. This value will 
decrease with rising T, vanishing around 30 mK. How- 
ever, from Eq. JHJ it can be seen that depends on T 
much more weakly (through the Boltzmann factors 
so that the IMT deficit can be reliably extracted even at 
temperatures where C cq is small or even zero. 

We stress that for quantum computing, it is the exis- 
tence of entangled eigenstates that matters, not the mea- 
sure of equilibrium entanglement [cf. the discussion a bit 
below Eq. since the system anyway must operate 
on a faster scale than the equilibration time. The cru- 
cial requirement is that the decoherence time T2 q of the 
entangled eigenstates exceeds the operation time. The 
presence of an IMT deficit shows that in our system T2 q 
is larger than min(fi,/A a / b , %j J) ~ 1 ns; its actual value 
can be determined using, e.g., the approach of |l4| . 

In conclusion, we have investigated two inductively 
coupled aluminum flux qubits with independently con- 
trolled bias fluxes. The impedance measurement tech- 



nique directly measures the low-frequency part of the 
system's magnetic susceptibility. The quantitative agree- 
ment between theory and experiment confirms that the 
system is in an equilibrium mixture of entangled states. 
Besides trying to scale the approach to larger qubit sys- 
tems, it will next be important to study the system's co- 
herence properties in detail, and to establish time-domain 
control enabling actual quantum gate operations. 
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